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HARMONIC ANALYSIS ON CAYLEY TREES II: THE BOSE 
EINSTEIN CONDENSATION 

FRANCESCO FIDALEO 



Abstract. We investigate the Bose— Einstein Condensation on non homoge- 
neous non amenable networks for the model deseribing arrays of Josephson 
junctions. The graphs under investigation are obtained by adding density zero 
5-H perturbations to the homogeneous Cayley Trees. The resulting topological 

%-H model, whose Hamiltonian is the pure hopping one given by the opposite of 

■^^ the adjacency operator, has also a mathematical interest in itself. The present 

paper is then the application to the Bose— Einstein Condensation phenomena, 
^L/ of the harmonic analysis aspects, previously investigated in a separate work, for 

Cm such non amenable graphs. Concerning the appearance of the Bose-Einstein 

Condensation, the results are surprisingly in accordance with the previous 

ones, despite the lack of amenability. The appearance of the hidden spectrum 
for low energies always implies that the critical density is finite for all the mod- 
els under consideration. First we show that, even when the critical density is 
r^ finite, if the adjacency operator of the graph is recurrent, it is impossible to 

exhibit temperature states which are locally normal (i.e. states for which the 
local particle density is finite) describing the condensation at all. A similar sit- 
r^ nation seems to occur in the transient cases for which it is impossible to exhibit 

I I locally normal states ui describing the Bose— Einstein Condensation at mean 

particle density p{(^) strictly greater than the critical density pc- Indeed, it is 
Cm shown that the transient cases admit locally normal states exhibiting Bose- 

^ Einstein Condensation phenomena. In order to construct such locally normal 

Cm temperature states by infinite volume limits of finite volume Gibbs states, a 

Cm careful choice of the the sequence of the chemical potentials should be done. 

'~fi For all such states, the condensate is essentially allocated on the base-point 

'-ft supporting the perturbation. This leads to p{u)) = pc- We prove that all such 

»Y~\ temperature states are Kubo— Martin-Schwinger states for a natural dynamics. 

/*^ The construction of such a dynamics, which is a delicate issue, is also done. 

(N 



1. Introduction 



The present paper is devoted to the phenomena of Bose-Einstein Condensation 
C^ (BEG for short) on non homogeneous networks obtained by adding density zero 

perturbations to homogeneous Cayley Trees. Even if the unperturbed model al- 
ready exhibit BEC (cf. [3] ) , the arising picture for the perturbed situation deserves 
attention as we have already explained in [6], and we are going to study in details 
in the sequel. The present paper is then the second part of the previous one [6], 
in which the harmonic analysis aspects of the model were intensively studied. The 
reader is also referred to [8] for some very interesting pivotal amenable models. We 
are going to investigate the physical application to the BEC of the very fascinating 
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picture which arises from the previous mathematical investigation. The physical 
underlying model concerns Bardeen-Cooper Boson pairs in networks describing ar- 
rays of Josephson junctions. The formal Hamiltonian is the quartic Bose-Hubbard 
Hamiltonian, given on a generic network G by 

Hbh = to ^ ni + ^ Aijiyn^fij - Joa\aj) . 

ieVG iJeVG 

Here, VG denotes the set of the vertices of the network G, al is the Bosonic creator, 
and rii — a^Oi the number operator on the site i € VG (cf. [4]). Finally, A is the 
adjacency operator whose matrix element Aij in the place ij is the number of the 
edges connecting the site i with the site j (in particular it is Hermitian). The reader 
is referred to the seminal paper [2] in which the theory of the superconductivity 
(called BCS Theory in honor of the authors) has been firstly discussed. It has been 
argued in [5] that, in the case when to, and V are negligible with respect to Jq, the 
hopping term dominates the physics of the system. After neglecting such terms, 
putting for the resulting effective coupling constant J, in principle different from Jq, 
J — 1, and finally passing to the one-particle space, the model under consideration 
becomes the so-called pure hopping ones described by the pure hopping one-particle 
Hamiltonian which assumes the form 

(1.1) H^\\A\\1-A. 

Here, A is simply the Adjacency of the fixed graph G, acting on the Hilbert space 
l?{VG). The constant \\A\\ is added just in order to have H > and does not affect 
the analysis at all. The reader is referred to [6] for the systematic investigation of the 
mathematical properties arising in non amenable cases considered here consisting 
of negligible additive perturbations of Cayley Trees. Several interesting amenable 
models are treated in details in [8]. The reader is also referred to [12] for some 
promising experiments for the pure hopping model on the Comb and Star Graphs 
(see Fig. 5) pointing out an enhanced current at low temperature which might be 
explained by condensation phenomena. 

One of the first attempts to investigate the BEC on non homogeneous amenable 
graphs, such as the Comb Graphs, was made in [5]. In that paper, it was pointed out 
the appearance of a hidden spectrum in the low energy part of the spectrum, which 
is responsible of the finiteness of the critical density. In addition, the behavior of the 
wave function of the ground state, describing the spatial density of the condensate, 
was also computed. Some spectral properties of the Comb and the Star graph were 
also investigated in [1] in connection with the various notions of independence in 
Quantum Probability. In that paper, it was noted the possible connection between 
such spectral properties and the BEC. 

The systematic investigation of the BEC for the pure hopping model on a wide 
class of amenable networks obtained by negligible perturbations of periodic graphs, 
has been started in [8]. The emerging results are quite surprising. First of all, 
the appearance of the hidden spectrum was proven for most of the graphs under 
consideration. This is due to the combination of two opposite phenomena arising 
from the perturbation. If the perturbation is sufficiently large (surprisingly, in some 
examples, it is enough that the perturbation is indeed finite), the norm \\Ap\\ of the 
Adjacency of the perturbed graph becomes larger than the analogous one ||A|| of 
the unperturbed Adjacency. On the other hand, as the perturbation is sufficiently 
small (i.e. zero-density), the part of the spectrum (T{Ap) in the segment {\\A\\, \\Ap\\] 
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does not contribute to the density of the states. This aUows to compute the critical 
density pdl^) at the inverse temperature /3 for the perturbed model by using the 
integrated density of the states F of the unperturbed one, 

^ ' P'^^P'- J e/^(a.+(||A^||-||A||)) _ 1 • 

The resulting effect of the perturbed model exhibiting the hidden spectrum (i.e. 
when ll^pll — ll^ll > 0) is that the critical density is always finite. The bridge 
between the Mathematics and Physics can be easily explained as follows. Consider 
the Bose-Gibbs occupation number (at the inverse temperature /3 > and chemical 
potential /i < 0) at low energies for the Hamiltonian (1.1). After using Taylor 
expansion, one hcuristically gets 

^,(H^,i)_j « [PiH ^/)]-i = [m\A\\ A*)/ - A)]-i ^ -^Ra{\\A\\ m) • 

Then the study of the BEC is reduced to the investigation of the spectral properties 
of the resolvent i?^(A), for A k, \\A\\, the latter being a very familiar object for 
mathematicians. Due to the particular form (1.1), as the BEC phenomena are 
connected with the spectral properties for small energies of the Hamiltonian, we 
can reduce the matter to the investigation of the spectral properties of the adjacency 
operator A of the network close to \\A\\. Another relevant fact connected with the 
introduction of the perturbation, and thus with the lack of homogeneity, is the 
possible change of the transience/recurrence character (cf. [11], Section 6) of the 
adjacency operator. It has to do with the possibility to construct locally normal 
states exhibiting BEC. As explained in [8], the last relevant fact is the investigation 
of the shape of wave function of the ground state of the model, describing the 
spatial distribution of the condensate on the network in the ground state of the 
Hamiltonian. From the mathematical viewpoint, this is nothing else the Perron 
Frobenius generalized eigenvector of the Adjacency (cf. [11]). It appears then clear 
that the physical and the mathematical aspects of the topological model based on 
the pure hopping Hamiltonian (1.1) are strongly related. 

Summarizing, the new and very surprising phenomena are the following. First 
the finiteness of the critical density is connected with the appearance of the hidden 
spectrum above the zero of the energy. It can be computed in a quite simple way 
for many amenable and non amenable models by using the secular equation, see 
Section 6 of [8] and Section 3 of [6] . This leads to the fact that the BEC can appear 
also in low dimensional cases rf < 3. In addition, not for all the models with finite 
critical density, it is possible to exhibit locally normal states (those for which the 
local particle density is finite in the infinite volume limit) describing BEC. Due to 
non homogeneity, particles condensate also in the configuration space, that is the 
system undergoes a " dimensional transition" . This phenomenon is connected with 
the transience /recurrence character of the adjacency operator of the network (the 
last describing the Hamiltonian of the system). Another new aspect is the difference 
between the geometrical dimension of the network and the growing of the wave 
function of the ground state. The last is called Perron-Frobenius dimension dpp 
for volume growing networks and is computed through the behavior of ^^-norms of 
the finite volume Perron Frobenius eigenvectors, all normalized to 1 at a common 
root. If the geometrical dimension da is greater than dpp, then it is impossible 
to exhibit locally normal states describing the BEC, whose mean density of the 
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particles is greater than the critical one. Surprisingly, there are non homogeneous 
networks for which pdP) < +00 and do — dpp (i.e. the comb graph N H Z^, 
see Fig. 1). In addition, there are also simple non homogeneous graphs exhibiting 
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Figure 1. The Comb Graphs N H Z'' and Z'' H Z. 

infinite critical density, but for which it is yet possible to exhibit locally normal 
states describing BEC. Such a very intriguing situation is summarized as follows.^ 
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A similar situation also happens for some non amenable networks. The networks 
under consideration in the present paper are density zero additive perturbations of 
exponentially growing graphs made of homogeneous Cayley Trees, see Fig. 2. The 
models treated in the present paper are the perturbations 0*^'', 2 < q < Q, and 
EII'3, of the homogeneous Cayley Tree G'^ along a subtree isomorphic to G"^, and N 
respectively, see below, see also [6] for the precise definition of such models and for 
further details. Quite surprisingly, the emerging mathematical situation described 
in [6] for such non amenable examples, is similar to that for the amenable graphs 
treated in [8]. Namely, we still find hidden spectrum, and the transience/recurrence 
character is determined by the transience/recurrence character of the negligible sub- 
graph supporting the perturbation. Finally, the infinite volume Perron -Frobenius 
eigenvector is exponentially decreasing far from the base-point supporting the per- 
turbation. Namely, concerning the wave-function of the ground state, we have 



A H B is the comb-shaped graph whose base-point is A, pc is the critical density, R/T denotes 
the recurrence/transience character of the Adjacency, BEC (p— BEC) denotes the existence of 
locally normal states exhibiting BEC (exhibiting BEC at any mean density p > pc). 
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precisely the same situation arising for the amenable situations treated in [8] (see 
Fig. 5). 

The first part of the present paper is devoted to describe general properties of 
(bounded) Hamiltonians on networks equipped with an exhaustion such that the 
density of the eigenvalues of the underlying Hamiltonian converges in the infinite 
volume limit to a cumulative function called Integrated Density of the States. It 
includes as particular cases those arising by negligible perturbations of Caylcy Trees 
as shown in [6]. Then we describe the general statistical properties of such models. 
In particular, a careful definition of the mean density and the possible amount of 
the BEC condensate should be done, because of the non negligible effects of the 
boundary terms due to the possible lack of amenability. We pass to describe the 
existence of the dynamics associated to the pure hopping Hamiltonian (1.1) for the 
perturbed Caylcy Trees. This is a delicate issue which is done in Section 4 for all the 
cases under consideration. After that, by using the results in [6], in the transient 
case we write down the formula (cf. (4.8)) for the natural class of states describing 
BEC (i.e. presenting a non negligible amount of condensate in the ground state 
wave-function). Such states are KMS (i.e. equilibrium) states for the dynamics 
previously investigated. As for the previous amenable models treated in [8], such 
states are infinite volume KMS states whose mean density coincides with the critical 
one. This simply means that the condensate is allocated on a subgraph whose 
volume growth is negligible with respect to that of the unperturbed one. We end 
the paper by investigating the infinite volume behavior of finite volume Gibbs states. 
Due to the lack of amenability, such an investigation presents some very delicate 
technical difficulties arising by the unavoidable boundary effects. Modulo such 
difficulties, the emerging situation seems to be in accordance with that described in 
[8], and summarized as follows. First, relatively to the general recurrent situation, 
it is indeed impossible to construct locally normal states exhibiting BEC at all, 
even if the critical density is finite. In addition, in the all transient cases under 
consideration, it seems to be impossible to exhibit locally normal states describing 
BEC at the mean density p > Pc{/3), (3 = 1/T. Finally, in the condensation regime 
where for the mean density p{uj) of all such states p{uj) — Pc{/3), a careful choice of 
the sequence of the finite volume chemical potentials /x(A„) — > should be done. 
By considering the geometry of the graph G and the £^-behavior of the ground 
state wave function ||f [a|| i such a choice corresponds to the local density of the 
condensate 

C.(A).D^^O, 

in the infinite volume limit A ^ G, see Proposition 2.3. This is nothing but the naive 
explanation of the fact that p{uj) = Pc{(3) for the locally normal states describing 
the condensation regime in all the situations under consideration. 

2. PERTURBED CAYLEY TREES 

A graph (called also a network) X — {VX, EX) is a collection VX of objects, 
called vertices, and a collection EX of unordered lines connecting vertices, called 
edges. Usually, the edges describe the physical interaction between vertices. In 
our model describing the pure hopping interaction, multiple edges as well as self- 
interactions, are also allowed. Let E^y be the collection of all the edges connecting 
X with y. We have E^y ~ Ey^. Let us denote hy A — [Axy]x,yex, x,y e VX, the 



adjacency matrix of X: 
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-^xy \^xy\ 



As explained in [6] , all the geometric properties of X are encoded in A. For example, 
X is connected if and only if A is irreducible. Setting 

deg := sup deg(x) , 

x£VX 

and 

Dxy := deg{x)5^^y . 

We have Vdeg < ||yl|| < deg, that is A is bounded if and only if X has uniformly 
bounded degree. For connected networks (or on each connected component) we can 
define the standard distance 

(2.1) d{x,y) := {mm £{TT{x,y)) \ 7r(a;, y) a path connecting x, y} , 

£{tt) being the length of the path tt (i.e. the number of the edges in tt). The 
adjacency operator considerably differs from the Laplacian A = A — D, whenever 
X is non homogeneous. On the other hand, it was proven in Theorem 7.6 of [8] 
that the value of the critical density does not change under negligible perturbations 
on a fixed homogeneous graph, if the Hamiltonian of the model is the (opposite of 
the) Laplacian. Thus, the main object for the investigation of the BEC for the pure 
hopping model of inhomogeneous networks will be the adjacency operator. In the 
present paper, all the graphs are connected, countable and with uniformly bounded 
degree. In addition, we deal only with bounded operators acting on £'^{VX) if it is 
not otherwise specified. 

Let B he a. closed operator acting on ("^{VX), and A e P(B) C C the resolvent 
set of i3. As usual, 

i?B(A):=(AI-S)-i 

denotes the Resolvent of B. 

Fix a bounded matrix with positive entries B acting on £'^{VX). Such an opera- 
tor is called positive preserving as it preserves the elements oi £'^ {V X) with positive 
entries. If B is selfadjoint, B is positive if {Bu, u) > for each u e £'^(VX). 

A sequence {v{x)}xevx is called a (generalized) Perron- Frobenius eigenvector 
if it has positive entries and 

J2 B^yviy) ^ \\B\Hx) , xeVX. 
yevx 

If such a vector is normalizablc (i.e. if it belongs to £'^{VX)) it is a standard 
i?^(V^X)-vector, otherwise it is a weight. In all the cases we simply call Pcrron- 
Frobenius eigenvector such a sequence by dropping the word generalized. 
Suppose for simplicity that B is selfadjoint. It is said to be recurrent if 

(2.2) lim {RbWSx,S^) =+oo. 

H\\B\\ 

Otherwise B is said to be transient. It is shown in [11], Section 6, that the re- 
currence/transience character of B does not depend on the base-point chosen in 
computing the limit in (2.2). The Perron-Frobenius eigenvector is unique up to a 
multiplicative constant, if X is finite or when B is recurrent, see e.g. [11]. In gen- 
eral, it is not unique, see the end of the present section for the case relative to the 
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Cayley Trees. We say that an operator B acting on 1^{VX) has finite propagation 
if there exists a constant r = r(B) > such that, for any x E X, the support of Bv 
is contained in the (closed) baU B{x,r) centered in x and with radius r. It is easy 
to show that for the adjacency operator if A on X, then A*' has propagation k for 
any integer fc > 0. 

Let X be an infinitely extended graph with an exhaustion {A„}„gN which is 
kept fixed during the analysis. The definition of the the integrated density of the 
states (cf. [10]) of a selfadjoint operator B E 'B{£'^(VX)) given in [7] works for all 
the situations for which it is meaningful, and then is useful for the purposes of the 
present work. Indeed, consider on 'B{i'^(VX)) the state 

Pn being the selfadjoint projection onto £^(yA„). Define for a bounded operator 
B, 

(2.3) T(B):-linir„(B), S G 2), , 

n 

where the domain 2),- is precisely the linear subspace of 'B{i'^{VX)) for which the 
limit in (2.3) exists. In addition, define for a bounded selfadjoint operator B, 

(2.4) T^ifiB)) := lim T„(/(P„i?P„)) , / G C(a(S)) , 

n 

provided such a limit exists. The domain D^-b C C* (B) is precisely the linear 
subspace of the unital C*-algebra C*{B) C 'B{i'^{VX)) generated by B, for which 
the limit in (2.4) exists. Notice that the definition of t, t^ might depend on the 
chosen exhaustion which we keep fixed during the analysis. 

Suppose now that T)^b = C*{B). Then it provides a state on C*{B) and, by 
the Riesz-Markov Theorem, a Borel probability measure /is on cr(-B). Thus, there 
exists a right continuous, increasing, positive function Fb satisfying 

Fb{x) — 0, X < min(T(_B) ; Fb{x) = 1 , x> maxcr(i?) , 

such that 

IJ.B{i—oo,x])=FB{x), xeR. 
The previous described cumulative function Fb is precisely the integrated density 
of the states (IDS for short) associated to B, provided it exists for the chosen ex- 
haustion. When the graph is amenable and the operator B is of finite propagation, 
the definition and some of the main facts relative to the IDS considerably simplify 
as the boundary effects play no role in the infinite volume limit, see Theorem 2.1 
of [7]. 

Consider the graph Y such that VY = VX, both equipped with exhaustions 
{Ar„}„gN, {Yn}nefi such that VYn ~ VXn, n G N. The graph y is a negligible or 
density zero perturbation of X if it differs from X by a number of edges for which 

\{e,yeEXAEY\xeVXn}\ 
™ \VX^\ 

where EXAEY denotes the symmetric difference. To simplify matters, we consider 
only perturbations involving edges, the more general case involving also vertices can 
be treated analogously, see [8]. Suppose that 1){t^^) = C*{Ax) that is Ax admits 
the IDS. 

An analogous result can be provided for the Comb Graphs in the transient situation. 
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Proposition 2.1. (Proposition 1.3 of \1]) Let Y he a negligible perturbation of the 
tree X. Then ©(t^^^) = C*{Ay), and 

T^nf{Ay))^T^Hf{Ax)). 

If |A| is sufficiently large, then it is possible to express the Resolvent of Y in 
terms of the Resolvent of X by the Krein formula. Indeed, fix X as the reference 
graph and define Ax '■= A, Ay :— A + D, where D is the perturbation, which 
eventually acts on 31(D). Put, for A G C, 

(2.5) S{X):=DPRa{X)\^j^, 



where P := P-T^jjy, is the selfadjoint projection onto 3^(1?) C ("^(VX). The Krein 
formula assumes the form 

(2.6) Ra^{X) = Ra. (A) + RA{X){lpmvx) - SiX))-'DPRA^ (A) . 

The resolvent formula (2.6) holds true for |A| > ||A|| + ||I?||, and extends to any 
simply connected subset containing the point at infinity of C U {oo} made of the A 
for which Ip^2(yx) ~ 'S'(A) is invertible. Conversely, suppose that 1p£2(^Yx) — S{\) 
is not invertible for some positive number A > H^jcll- Then the norm ||^y|| of the 
perturbed graph Y is greater than that ||Ax|| of the unperturbed one X. It can be 
computed as 

||yly|| = sup{A > I Ip£2(yx-) — 5'(A) is not invertible} . 

Fix now our set-up which is a graph X equipped with an exhaustion {Xn}neN 
such that Ax admits the IDS. Let y be a negligible additive perturbation of X as 
explained above. First of all, note that Proposition 2.1 assures that Ay admits the 
IDS w.r.t. the exhaustion {KnlneN- In addition, we obviously have ||j4^|| > ||A^||. 
Define 

(2.7) S:=\\Ax\\-\\Ay\\. 

Denote Fx ■= -^||Ax||i-^x' ^y '■— ^\\Ay\\i-Ay^ where Nb is the IDS of the operator 
B. 

Proposition 2.2. (Corollary 2.6 of [Q\) We have 

Fy{x)^Fx{x + 5). 

Proof. The IDS of ||^jf j|I - ^x and ||yly||I - Ay exist by hypothesis and by 
Proposition 2.1, respectively. We get, again by Proposition 2.1, 

j f{x)dFy(x) = T^^{f{\\Ay\\l-Ay)) = T^- {f {\\Ay\\l - Ax)) 

=T^Hf{\\AxU-Ax-5t))= f f{x-S)dFx{x)^ f f(x)dFx{x + S) . 



This leads to the assertion. D 

Suppose that \\A^\\ > ||A'^||. Then 5 < in (2.7). It has a very precise physical 
meaning as an effective chemical potential (cf. (1.2)), see Proposition 7.1 of [8]. 
Thus, the part of the spectrum (||Ay|| +S, ||Ay||] does not contribute to the IDS of 
Ay. Equivalently, the perturbed pure hopping Hamiltonian ||Ay||I — Ay presents 
the hidden spectrum, see Proposition 3.2. The appearance of the hidden spectrum is 
then the combination of two different effects: the perturbation should be of density 
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zero in order not to affect the IDS (Proposition 2.1), but it should be sufficiently 
large in order to increase the norm of the perturbed Adjacency. 

Our framework will be the investigation of the occurrence of the BEG for per- 
turbed Cayley Trees, whose mathematical aspects are investigated in detail in [6]. 
Namely, we consider the graphs G'^^'' with q > 2, and q < Q < Q{q) where 



(2.8) 



Qiq) 



2v/g^^+i + /V9^^ 



1 



{[x] stands for the integer part of a positive number x), and H'^ with 2 < Q < 7. 
Such networks consist of the Cayley Tree G'^ of order Q, suitably perturbed along 
subtrees as explained in Fig. 2. One can prove (cf. Lemma 4.1 of [3], together 






Figure 2. The networks G^'^ H^ and G'^^^^ 



with Proposition 2.2 of [6]) that the Adjacency of G*^ admits the IDS w.r.t. the 
standard exhaustion made of all the balls of radius n centered in a fixed root. 
In the above notations, the Laplace transform <I> = L[Fqq] (or equivalently the 
one-particle partition function in physical language) of Fqq is given by 



(2.9) 



^ ^ fe=ln=l 



■k^-4PVQ^sin- 2(ffT) 



This, together with Proposition 2.2, leads to Fy{x) = Fqq(^x + ||^gq|| ~ ll^i'll)! 
where Y is G'5''' or H'^. As \\Ay\\ — II^g'^II > Oj they all present hidden spectrum. 
We also have that the graphs G'^'' and H'^ are all transient except G'^'^^ see [6]. 

We briefly report the description of the Perron Frobenius eigenvectors v given 
in [6] for the graphs under consideration. Let S be the base-point of the per- 
turbed graph Y under consideration, supporting the perturbation. The exhaustion 
{A„}„gN has an obvious meaning by considering the part Sn of 5* in B^ , the last 
being the ball of radius n of G*^, centered in a root 0, kept fixed during the analysis. 
It was shown in [6] that 



(2.10) 



v{x) = a(\\AY\\f^'''^'>w{y{x)) , xeVX. 



Here, a(A) is the function given (4.2), d is the distance given in (2.1), y{x) is the 
unique point on S such that d{x, S) ~ d{x, y(x)) (cf. Lemma 4.2 of [6]), and finally 
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w is the Perron- Frobcnius cigenveetor of PRa .q (ll^r ||)^ {P ^ Pp{s)) given by 
(2.11) w{x) = {l-a{\\AaQ\\))d{x,0) + l, xeS^N, 



w{x)=[l + - rf(x,0))(g-l) — ^, xeS* 

described in Theorems 5.2 and 6.3 of [6] for M9 and G''''', respectively. Concerning 
the behavior of the ^^-norms of the restrictions of u, we have the following 

Proposition 2.3. For the Perron-Frobenius eigenvector v described by (2.10) and 
(2.11), we get 

(2.12) ]imtki^^J^ = . 

|A„| 

Proof. First of all we note that, in all the situations, 

yes 
This leads to 

II II 2 II II 2 

pr^2(A„)|| \Sn\ , II ,,, |nMK2|Pr<?2(S'„)|| \Sn\ 



l,^^<(l + ||«.„(p,-||))||) „,_, 



By taking into account the form (2.11) of w, we easily obtain the following estimates. 

|2 
|2, 



Concerning H'5 we have w[^2(5^) « n^, which leads to 



rr^=(A„)|| I'S'nl ^ n^ 



|A„| 0" 

For the case G'^''?, q = 2 leads to w{x) — 1, identically on S. Concerning the case 
g > 2, consider the function 

It is not difficult to see that /(I) < 1 and /'{£_) < 0, ^ > 1. This leads to w{x) < 1 
whenever x G <S'\{0}. We finally get 



<(!+ R.^i\\Ay\\))rL^^^^O. 



|A„| -- ' "-"-«^" ""^^11^ \B^\ - g- 

D 

We end the present section by noticing that it is easily seen that the Perron- 
Frobenius eigenvector is no longer unique, in general for infinitely extended net- 
works, and in particular for the situation considered here. In fact, on G*^, Q > 2, 
'Pq,xo{^)j ^ £ C^ are Perron-Frobenius eigenvectors for the Adjacency, for any 
fixed root xq E G'^. The same happens for G'''^, q > 2 where in (2.10) it is enough 
to choose w{x) — ipg^^^^lx), xq £ S ^ G"^ is any fixed root, where 

(Pq.y{x)= U + - d{x,y)] (q-l)-^^^ , a;,yeG«. 

The choice of the normalization at 1 on a fixed root (which is denoted by in our 
framework) selects in a unique way the Perron-Frobenius vector among the class 
previously described. For the networks G^ = Z and G'^'^, the Perron-Frobenius 
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eigenvector for the Adjacency is unique (up to a multiplicative constant) thanks to 
the fact that the Adjacency is recurrent. It is expected that it is unique also for 
H'5 , thanks to the fact that 5 ^ N has a unique end at infinity. 

3. STATISTICAL MECHANICS ON INFINITELY EXTENDED NETWORKS 

In order to investigate the statistical properties of the pure hopping model de- 
scribing Bardeen Cooper pairs on arrays of Josephson junctions on non homoge- 
neous networks we report some standard notions useful in the sequel. 

Let (2t, a) be a dynamical system consisting of a non Abelian C*-algebra and 
a one parameter group of ^-automorphism a.^ The state tp on the C*-algebra 
*8 satisfies the Kubo-Martin-Schwinger (KMS for short) boundary condition at 
inverse temperature /3 e IR\{0} w.r.t the group of automorphisms {rtj^gji if 

(i) for every A, _B e *B, i h^ (p{Aat{B)), 1 1-^ (p{at{A)B) are continuous; 
(ii) for each f el), 

(3.1) J ^{AaAB))f{t) dt^ I ^(at{B)A)f{t + ip) di , 

where " ^ " stands for the Fourier transform and T> is the space of the 
infinitely often differentiable, compactly supported functions on M. 

The following facts are well known. First, a KMS is automatically invariant w.r.t. 
the automorphism group at- Second, the cyclic vector Q.^p of the Gelfand-Naimark- 
Segal (GNS for short) quadruple (7r^,?{^, C/y,r2^) is also separating for 7r^(S)". 
Denote a'^ its modular group. According to the definition of KMS boundary con- 
dition, we have 

af OTT^ — TT^ O a-i3t ■ 

We refer the reader to [4], Section 5.3 for equivalent formulations of the KMS 
condition, proofs and further details. 

The C* -algebras considered here are those arising from the Canonical Commu- 
tation Relations (CCR for short). Namely, let f) be a pre-Hilbert space equipped 
with the non degenerate inner product (•, •). Denote 1} its completion. Consider 
the following (formal) relations, 

(3.2) aif)aHg)-aHg)aif) = {gJ) f,gei). 

It is well known that the relations (3.2) cannot be realized by bounded operators. 
A standard way to realize them is to look at the symmetric Fock space 9^+(f)) on 
which the annihilators a(/) and creators a^f) naturally act as unbounded closed, 
adjoint each other (i.e. a(f)* = a^f)) operators. This concrete representation of 
the CCR is called the Fock representation. An equivalent description for the CCR 
is to put on 9^+ (5), 



(3.3, *(/, :. 'An±^ 

and define the Weyl operators W{f) :— exp «$(/). The Weyl operators are unitary 
and satisfy the rule 

(3.4) W{f)W{g)=e''''^W{f + g), .f,gei). 



No a-priori regularity assumption is made on tlie automorpliisms group {at}t 



6K 
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It is well known that the representation on 3^+(i)) of 2IT(f)) described above, is faith- 
ful. We adopt the Weyl algebra as the definition of the CCR algebra, the last being 
the abstract C*-algebra 2U([}) generated by the operators {Vr(/)}/eij satisfying 
the relations (3.4). Sometimes, it is also useful to compute directly the expectation 
values of combinations of products of annihilators and creators (equivalently field 
operators (3.3)) in a given state. 

Let H he a self adjoint operator acting on i). Suppose that e**^f) C f). Then 
the one parameter group of automorphisms Ttf :— e**^/ defines a one parameter 
group of ^-automorphisms at of 2IT(f)) by putting at{W{f)) := W{Ttf). The latter 
is called the one parameter group of Bogoliubov automorphisms generated by Tt . 

A representation tt of the Weyl algebra 2U([}) is regular if the unitary group 
A € M I— )■ 7r(VF(A/)) is continuous in the strong operator topology, for any / G t). 
A state (p on 2U(f)) is regular if the associated GNS representation is regular. 

The quasi-free states on the Weyl algebra are those of interest for our purposes. 
Such states ui are uniquely determined by the two-point functions uj(a^ {f)a{g)), 
f,g G i). A standard textbook for CCR is [4] (cf. Section 5.2) to which the reader 
is referred for proofs, literature and further details. 

Let G be any graph. To simplify notations, we denote the sums on the vertices 
VS C VG of subgraphs S C G, and of VG itself, directly by sums on S or G, 
respectively. Consider a subspace f) oi £'^{G), which contains the indicator functions 
{6x I X G G}. A representation tt of the Weyl algebra 2IT(f)) is said to be locally 
norm,al (w.r.t. the Fock representation) H tt\w{P{A)) is quasi equivalent to the Fock 
representation of 2n(i?^(A)). A state on 2IT([)) is said to be locally normal if the 
associated GNS representation is locally normal. A locally normal state cp does 
have finite local density 

I I jeA 

even if the mean density might be infinite (e.g. lirasupf^^Q p\{ip) = +oo). Let 
A„ t G be an exhaustion of G, that is a sequence of finite regions invading the 
graph G, together with a sequence of states {oja„} on 2n(^^(A„)). 

Proposition 3.1. (Lemma 3.2 o/ [8]J Suppose that 

\iT[VLij{^^{a* {5j)a{5j)) = +oo 

for some j G G. Then uj(W{v)) := liTanLU\^{W{v)) does not define any locally 
normal state on CCR{i}), where () C l'^{G) is any subspace containing the finite 
supported sequences. 

Our purposes concerns the pure hopping model on the networks G*^'"^ and H'^ 
or, more generally on a uniformly bounded connected network X equipped with 
an exhaustion {X„} for which the Adjacency admits the IDS (cf. [7]). As previ- 
ously explained, all the physical properties of the emerging quasi-free multi-particle 
model are encoded in the one-particle pure hopping Hamiltonian. Namely, the 
graph X will be one of such networks, Y a density zero additive perturbation of X 
as described above, equipped with exhaustions {Xn\n&ii, {i^njnGN- We denote by 
G one of such graphs, equipped with the finite volume exhaustion {A„}„gN- The 
pure hopping Hamiltonian, after normalizing such that the bottom of the energy is 
zero, assumes the form 

nAaW-Ac- 
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The finite volume Haniiltonians will be _ff„ := PnHPn acting on ^^(A„), P„ = 
Pp{A„) being the selfadjoint projection onto £^(A„). 

Fix a general bounded positive Hamiltonian H E 'B{£'^{G)) admitting the IDS 
for a fixed exhaustion as for the pure hopping one.^ Consider the finite volume 
densities of the states Nn relative to iJ„, and that N relative to H. Define 

Eo(H) ■- lim f infsupp(A^„) ) =0 
A„tG y ^ J 

Em(H) :=inf supp lim Nn I = inf supp (A'^) , 

\A„tG / 

where the limit in the first equation exists by Lemma 3.4 of [8], and that in the 
second one is meaningful directly by the definition of the IDS. We always get 
Em{H) > Eo{H) = 0. 

Definition 3.2. If Eq{H) < Em{H) we say that there is a (low energy) hidden 
spectrum, see e.g. [5]. 

The following assertions can be shown for all the situations under consideration. 
First, H is continuous in (see also [8], Proposition 3.9) in all the situation con- 
sidered here. Second, in the unperturbed situation the hidden spectrum does not 
occur (by direct inspection of (2.9) for 0*5, or by Theorem 5.2 of [8] in the case 
of periodic amenable graphs). If Y is any perturbation of X, we have for the 5 
appearing in (2.7), 

E^{Hy) - ErniHy) = -Em(H) = \\Ax\\ - \\Ay\\ - <5 . 

The opposite of Em{H) then plays the role of a chemical potential in determining 
the critical density of the perturbed graphs under consideration. Indeed, consider 
the formula for the mean density at the inverse temperature /? of the particles (cf. 
[4]) for the model with Hamiltonian H > 0, and chemical potential /x < 0, 

f dN„(h) 

and the corresponding critical density 

,«,ffl:=/^.p„(AO). 
By taking into account Proposition 2.2, we get 

This means p^^ {/3) < p^^{f3). In addition, if there is hidden spectrum, p^^ (13) is 
always finite and p^^ {/3) < p^^{f3). It can happen that, after perturbing a graph 
whose critical density is infinite, we get finite critical density. As for the networks 
under consideration in the present paper, the unperturbed critical density is always 
finite (cf. [3]), the only effect of the appearance of the hidden spectrum on the 
critical density is that it decreases w.r.t. the unperturbed one. 

It is well known that the condensation regime is described by /i = 0, see e.g. 
Section 5.2 of [4], or Section 3 of [8]. Indeed, one starts from the Gibbs grand- 
canonical ensemble at finite volume A C G and fixed density. One fixes an inverse 



As usual wc normalize H such that G <^{H) is the bottom of the spectrum cr{H) of the 
Hamiltonian H . 
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temperature /3 > together with the particle density p, and determines the finite 
volume chemical potential /i(A) by solving 

(3.5) phA(3,K^))^P, 

see e.g. Fig. 3. The infinite volume limit can be investigated by means of the 




Figure 3 

reference exhaustion {A„}„gN. To take into account also the very different situation 
appearing in non homogeneous situation, we start by fixing any general sequence 
of chemical potential 

fi„ < \\Ag\\ - \\A^J ^ EoiH^J , 

which we can suppose to converge (eventually passing to a subsequence) to some /i. 
In the case when such a sequence is recovered by using (3.5), we put (!„ :— /i(A„). 
Since Eq{H{^^) -^ Eq{H) = 0, we get ^ < 0. The finite volume state with chemical 
density p at /3 is described by the two point function 

(3.6) oj^^iaHOaiv)) - ((e^(«-"-^"^) - l)-'^,v), 

where ^,?7 € ^^(A„). Thermodynamical states are then described as limits of the 
finite volume states as above. Concerning the infinite volume limit of the finite 
volume density, we report a result (cf. [5]), whose proof relics upon the analogous 
one, Proposition 3.8 of [8]. It describes in a rigorous way how it is possible to define 
the amount of the condensate in the infinite volume limit for an arbitrary sequence 
of finite volume chemical potential /x„ — > 0. Due to the general character of such 
a result, we suppose that _ff is a bounded positive Hamiltonian on a network G 
equipped with an exhaustion A„ such that H admits the IDS N (in the situation 



under consideration H 



\Ar 



Ar 



Denote as usual P„ 



P,2 



(A„), 



and 



Hn '■= PnHPn = -ffA„ • To avoid the almost trivial situation /i < 0, we restrict the 
matter to the condensation regime fi — 0. 

Proposition 3.3. Under the conditions listed above, the following assertions hold 
true for each sequence ii„ < Eq{H„) satisfying ii„ — > 0. 

It is also customary to fix tlic activity z := e'^'^, instead of the chemical potential. 
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(i) Suppose that H admits low energy hidden spectrum. Then for each contin- 
uous mollifier F which is identically in a neighborhood of 0, and 1 in a 
neighborhood of supp A^ , we have 



lira 

A„tG 



gP(X-fJ. 



g^diVH. (x) = pf(;3)<+cx3. 



(ii) Suppose that linix-ia N(x) = 0. Consider a sequence {F^}g^() of increasing 
continuous mollifiers, all vanishing in a neighborhood of and converging 
almost everywhere (w.r.t the measure determined by N ) to 1. Then 

lim lim / ^^l""] ANh,^ [x) ^ pf (/3) . 

In addition, if p^ {(3) < +oo, one can replace for the continuous mollifiers 
{-Fe}e>0; the condition of monotony with uniform boundedness. 



Proof. We start by noticing that 



Fix) 



and all the 



F.jx) 



are continuous in a com- 



mon neighborhood of all the supports of the measures determined by the cumulative 
functions Nh^ and Nh- 

(i) By the existence of the hidden spectrum, we have 

F{x) ,,, , , [ ANh{x) 



< 



g/3(a;-/^„) _ I 
1 



dNH„{x) 
1 



gf3(x-fj.„) _ I 



oPx _ 2 



efSx _ I 
F{x)dNHjx] 



^f)x _ ^ 



dNuix) 



since 



.,3(x-K„)_l 



,fi^_^ , uniformly on the support of F, and the second summand 
goes to zero because / f{x) dNn^ (x) -^ J f{x) ANh{x) for each continuous function 
/ in a common neighborhood of all the supports of the involved measures {d-/V„}„gN 
and AN. 

(ii) The same computation as before leads for each e, 



Feix) 



g^(a;-/j„) _ I 



ANhA^) 






dNuix) 



The proofs follows either by the Monotone or the Dominated Convergence Theorem, 

asO = iVH(0)=lim,^oA^H(2;). □ 

Thanks to the above result, the quantity 

Un := lim lim / —tt-, -. dNtj. (h) 

is well-defined and independent of the particular choice of the mollifiers as above. 
Such a quantity describes the amount of the condensate in the ground state. Indeed, 
if the sequence of chemical potentials are obtained through (3.5) by fixing the mean 
density p, we have 

p = no + p{f3,p,). 

Remark 3.4. 

(i) yU < if and only if p < pdfi)- 
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(m) For any /3 > 0, /i < 0, the sequence (3.6) converges pointwise to a state w, 
whose two-point function is given by 

Moreover, the density p{lu) of the state lu, defined by 



1 

I 



p{uj) :=hnisup--^ V WA„(a^((5j)a((5j)) 



satisfies p(w) — p{(3,p). 

{Hi) The transience of A is a necessary and sufficient condition (see Proposition 
4.4 and Remark 5.1) for the existence of locaUy normal states on 2n(f)) describing 
condensation phenomena, that is when for p, — Q^ for the pure hopping modcL 
{vi) Observe that ng > if and only if p > pdP)- Indeed, hq = whenever /i < 0, 
hence hq > ^ p = 0. As a consequence, < no = p — p{l3,p) ~ p — pc{l3), i.e. p > 
Pc{P). Conversely, p > Pc{P) implies /i = 0, hence no = p — p{l3, p) — p — pdfi) > 0. 
(v) Observe that no > can be obtained only if /i„ — > 0. Indeed, /^ = is a 
necessary condition for the occurrence of BEG. 

(vi) Observe that no > can occur also when the infinite volume limit of the 
two-point function leads to non locally normal states i.e. uja„{<i''' {Sj)a{Sj)) — > +oo. 
(vii) Observe that it is possible to have locally normal states exhibiting a non trivial 
amount of condensate even if p = Pc{/3), see Section 10 of [8], or Sections 4, 5 below. 

To simplify the notations, from now on we put l3 = 1 throughout the paper if it 
is not otherwise specified. We also put p{p) := p^{l,p), pa{p) :— p^'^{l,p), and 
finally pc := P'^(l, 0), H being the Hamiltonian of the model. 

4. THE EXISTENCE OF THE DYNAMICS AND THE KMS STATES EXHIBITING BEG 

Let Y be any of the perturbed graphs under investigation. In order to study 
states describing local properties, and in particular those which are locally normal, 
we start to consider a subspace f) C l?{Y) which contains all the canonical elements 
(5j, and which is stable under the dynamics generated by if := \\Ay\\ — Ay- To this 
end, we define 

(4.1) f) :=span {e**-f^(5j I teM,.? e r} 

for the algebraic span (no closure in P{Y) = l'^{Gi^)) of all the elements of the form 
e^^^Sj. By construction e**^f) C f). For our aims, we need also that the pairing 
/i G f) n- (/i, v) is well defined, when v is the Perron Frobenius vector, described in 
[6], for the adjacency matrix Ay in each of the cases under consideration, see (2.10) 
and (2.11). For u G i^{Y) we also define |-u| € £'^{Y) by putting |u|(a;) := \u{x)\, 
x€Y. 

Lemma 4.1. Let f :~ X^nS'^"-^" ^^ ^ function analytic in a neighborhood of 
(T{Ay) such that a„ > 0, n G N. Then f{Ay)v = f{\\Ay\\)v. 

Proof. As the adjacency matrix has positive entries, all the involved series have 
positive entries. Thus, by the Monotone Gonvergence Theorem we get 

+ 00 +00 /+00 \ 

f{Ay)v = Y, anA^v = J2 anWAyirv = ^ aJAy]]" V = f{\\Ay\\)v . 

n=0 n=0 \n=0 ) 

U 
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The following functions, defined in R\[-2y/Q - l,2y/Q - 1], 



(4.2) a(A) := 



1-a/I- 



4(Q-i) 
A2 



2(Q-1) 
A 



(4.3) fi{\) := 



'-2 + qVi-^^ 

2(Q-i) 
A 



are useful in the sequel. If 5 C G*^, it was shown in [6] that, for the situation under 
consideration, the operator S{X) in (2.5) describing the Resolvent of Ay through 
the Krein formula (2.6) assumes the form 

eventually acting on i^{S).^ In the cases under consideration when Y is G'^''^, H'5, 
5 ^ G^ (where G^ ~ Z) and S' ~ N, respectively. In all such situations, (2.5) is 
proportional to the convolution by (4.2), and for z £ C\[—2^/Q— 1, 2^/Q — 1] it 
assumes the form 

(4.4) iS{z)u) (x) = ^^"^ '. ^ , 

provided at least that the series in (4.4) converges pointwise. Here, d is the dis- 
tance function given in (2.1), and a{z) and /x(z) are the analytic continuations 
in C\[— 2^/Q — 1, 2y/Q — 1] of the corresponding objects given in (4.2) and (4.3), 
respectively. It was shown in [6] that (4.4) is meaningful and defines an element 
of £^(G'5), at least when A > A* = ||Ay||. Now we show that the series in (4.4) 
converges and defines an element in ^^(G*^) for z e C such that l^l is sufficiently 
large. 

Lemma 4.2. Let a{z) and ii{z) as above, and r > |l^y|l > 2^/Q — 1. Then 

\a{z)\ a{r) 



r 



{z\\z\=r} \^i{z)\ fl{r) 



Proof. Put w := 2^Q — 1/ z. Disregarding the constants and the quantities whose 
modulus is fixed, we have to maximize and minimize on circles of fixed radii less 
than 1, the functions 

2 2 

F{w):= 1- Vl-w^ , G{w):^ h+y/l-w^ 

where w is in the open punctured disk int(Di)\{0} of radius 1 centered in the origin 
of the complex plane, and b = (T^^ see Fig. 4. First of all we note that, for the 

analytic continuation Vl — w'^ of Vl — x^ in the cut plane C\((— oo, — 1]U[1, +oo)), 
we have 



(4.5) Vl - i-w? = Vl - ^2 , Vl - {^? = VT^ 



With an abuse of notations, we are dropping the dependence on the subset S on which is 
allocated the perturbation. 
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Figure 4 

After some tedious computations, the functions F, G are expressed by functions /, 
g on the unit circle, given by 

/(^) := 1 + hi-d) - V2^h{d) + 1 - p2 cos 2d , 

g{'&) := b^ + hi-d) + by/2^Jh{d) + I - p^ cos 27? . 

Here 

hi-d) = v/(l + /52)2 - V C0S2 d , 

< yO < 1 is fixed, and by (4.5), we can reduce the matter to < ^ < tt/2. After 
some straightforward computations, we get that /'(i?) and .g'(z?) are never zero for 
< I? < tt/2 so the the functions attain the extrema either in or in tt/2. We get 
for /, 

v/7(0) = i-v/i-p'>v/i + p'-i = yiW2), 

because 4 > 2(1 + -y/l — p^). Concerning g we have, 

VgW) = b+^l-p^<b+^l + p^ = ^g{TT/2) . 
Collecting together, we get with p = 2^/Q — 1/r, 

P ^ 1 ^^^ IJW < niaxo<^<^/2 y/m ^ l_ //(O) ^ aM < p 
Q 0<™<V2 V 5('?) " mino<^<^/2 VsM Q^ 9{(^) K^) ~ 
and the proof follows. D 

Denote D,, C C the closed disk centered in the origin of radius r. 

Proposition 4.3. There exists r > \\Ay\\ > 2\JQ — 1 such that for S{z) in (4.4), 
we get ||S'(2;)|| < 1 for each z e C\Dr- Thus, for z e C\Dr, the resolvent of the 
Adjacency of Y can he written by the Neumann series 

(4.6) RaA^^) - Ra^^z) + Ra^^Jz) [Y^SizA P,2(s)Ra,Jz) . 

\k=0 / 
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Proof. Let u ~ ui — U2 + »(u3 — U4) be the decomposition of m G £"^{3) in positive 
elements. One has \\uj\\ < \\u\\, j ~ 1,...,4. By taking into account (4.4) and 
Lemma 4.2, we get 

4 

{S{z)u,S{z)u) < J2{\S{z)ur\,\S{z)u,\} 

4 

<Y,{S{\z\)u,,\S{\z\)u,)<W\\S{\z\)r\\ur 

As ||S'(A)|| is continuous and decreasing (cf. Lemma 3.1 in [6]) with hmA-f+co ll-^CA)!! = 
0, there exists r > \\Ay\\ such that |j5'(r)|j < 1/5. CoUecting together, we have in 
the complement of D,. , 

ll^Wll<4!15(r)|l<4/5<l. 

In any simply connected open subset of C U {oo} containing the point at infinity 
for which l£2(-5) — S{X) is invertible, Ray{z) assumes the form 

RA^iz) = Ra„q (z) + Ra,.q {z)Rs(z){'^)Pi2i^s)RA^,Q {z) , 
and the proof follows from the first half, by using the Neumann expansion of 

Rs[z){l)- □ 

Consider the function 

l^ e^ — 1 X ' 

It is bounded on [0, +oo). In addition, ^J-_^ = f(x) + i. It provides the comparison 
between the resolvent of A and the functional calculus of (e^ — I)~^ associated to 
the Bose-Gibbs occupation number for the pure hopping Hamiltonian. 

Proposition 4.4. Under the above notations, the following assertions hold true. 
(i) Ay is transient if and only if 1} d D ((e-^ — I)^-^/^) . 
(ii) For h (^ i:} we have X^ksf \^i^)\'^i^) < +oo. 

Proof. It is enough to prove the assertions for the generators h = e^^^5j of (). 

ii) By using the splitting described by (4.7) we get that w is in I> ((e^ — I)^^/^) 
if and only if w is in D(i?yi(|| A||)^/^). Let j be any element of Y . Obviously, 5j e [} 
and, by hypothesis, 5^ G T) {{e" - t)-^/"^). So 5, e 1){RA{\\A\\y/^) and Ay is 
transient. For the reverse implication, we compute for A > ||Ay||, 

(i?^, (A)e^*^,5,, 6**^,5,) =(i?^,(A)e^*(ll^-ll^-^-),5„e^*(ll^-ll^-'4-)(5,) 

= (e^*(ll-^-ll^-'^-)i?A,(A),5„e^*(ll^-ll^-'^-)(5,) 

= {Ray{X)Sj,6,) t (RA^iWAylDS^^Sj) 

which is finite if Ay is transient. Thus, e'*"Sj € T> ((e^ - I)"!/^). 

{ii) Fix R > r, where r is the number appearing in Proposition (4.4), together 
with the circle Cr C C of radius R centered in the origin. By using the Neumann 
expansion of Rs(^z){^) on C/j and the fact that Ra q (z) itself is the convolution of 
the function a{z)^ / ii{z) (cf. (7.6) in [9]), we show that the corresponding resolvent 
Ray {z) is a combinations of convolutions and infinite (absolutely converging) sums 
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of powers of the convolution by the function a{z)^ / ^{z), see (4.6) and (4.4). Let 
Cn be counterclockwise oriented. Thus, thanks to Proposition (4.4) as explained 
above, together with Lemma 4.1, we get 

(|e^*%|,^')=( ^£ e^*("-^-"»-)i?^,(z)<5, dA ,«^ 

<R{\RaAz)5j\,v) <R{Ra^{R)S,,v) 

R{Sj,v) Rv{j) 



--R{Sj,Ra^{R)v) 



R-\\Ay\\ R-\\A^ 



D 



Now we pass to describe a class of locally normal states exhibiting BEC, which 
are indeed KMS for the natural dynamics generated on 2n(f)), f) given in (4.1), 
by the Bogoliubov automorphisms whose generator is the (second quantization of) 
the pure hopping Hamiltonian. To this end, we consider the dynamical system 
(2n([)),at), where 22J(f)) is the C*-algebra made of the Weyl CCR algebra on (), 
and at is generated by the Bogoliubov automorphisms TfU := e**^u, 

atiWiu))=W{Ttu), ue\), 

which are well defined as e**^f) = f) by construction. 

By taking into account Proposition 3.1, the necessary condition for the existence 
of locally normal states exhibiting BEC is the transience of Ay- For the pure 
hopping model, the transience of the Adjacency is indeed also sufhcicnt as we are 
going to see. For the rest of the present section we limit our analysis when Y 
is either U^ , 3 < Q < 7, or G^'? where g > 3, and g < Q < Q{q) with Q{q) 
given in (2.8), see Theorems 4.1, 5.3, 6.5 of [6], together with Proposition 3.1. Fix 
D > and define states on the Weyl CCR algebra 2II(()) uniquely determined by 
the two-point function 

(4.8) ujD{a\ui)a{u2)) -.^ {{e" -iy^ui,U2) + D{ui,v){v,U2} , ui,U2el)- 

The states ujd are well-defined (cf. [4], Section 5.2) because of Proposition 3.1. In 
the case of quasi-free states on 2n(i}), it is also customary to provide the description 
directly in terms of their values on the Weyl unitaries. We get (cf. Section 5.2.5 of 

[4]), 

ujd[W[u)) ^ e 4 e 2 ^ ^ e f) . 

Note that for quasi-free states, and then in our situation, the KMS condition can 
be checked directly in terms of ij£){a^ {ui)a{u2)) , at least when 

1 1-> ujD{a\ui)a{TtU2)) 

is bounded and continuous. The boundedness automatically follows by the defini- 
tion of the two-point function thanks to transience. The continuity can be easily 
checked by Lebesgue Dominated Convergence Theorem by taking into account again 
Proposition 3.1. 

Theorem 4.5. For each D > 0, the states ljd given in (4.8) are KMS (at inverse 
temperature j3 = 1) for the dynamics generated by the Bogoliubov automorphisms 
{e'"* I teM}. 
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Proof. It is easily seen by Proposition 3.1 that ujd{o.' {ui)a{TtU2)), wi3(a(TtU2)a' (ui)) 
are bounded and continuous (in t) . The proof follows by an elementary change of 
variables in (3.1) as / e D. D 

The second addendum in (4.8) describes the amount of particles which conden- 
sates in the ground state whose wave function is nothing but the Perron-Frobenius 
weight V. In the finite region A, such a density of the condensate is roughly given 

by 

|2 



(4.9) 



Cd{A) 



1 

W\ 



J2d{S^,v}{v,6,)^D- 



xeA 



|A| 



Even if the model exhibits a non negligible amount of the condensate, C£)(A) — >■ 
as A t i^, see Proposition 2.3. Formula (4.9) explains what happens in the non 
homogeneous pure hopping model. Particles condensate even in the configuration 
space. The amount of particles that the system can accommodate in the subnetwork 
made of the support of the perturbation is governed by the ratio between the volume 
growth of the original network and the ^^-growth of the wave function of the ground 
state. In all the situations into consideration in the present paper, Cjj{A) goes to 
as A approaches to Y. This is due to the quite surprising fact that the Perron- 
Frobenius weight decreases exponentially far from the perturbed zone (cf. Fig. 5), 
see [6], or also [8] for several amenable cases. This is the naive explanation that the 




Figure 5. The wave function of the ground state for the Comb 
and the Star Graphs. 

mean density of the particles associated to all the states lod > exhibiting BEC is 
the critical one pc- 



5. INFINITE VOLUME LIMITS OF GIBBS STATES 

We start with a simple general result explaining that, for the pure hopping 
model, the transience is a necessary and sufficient condition for the existence of 
locally normal states exhibiting BEC if one disregards the dynamics. 

Remark 5.1. Let G be a uniformly bounded degree connected network. Suppose 
that its Adjacency A is transient. Let ()o := span {Sx \ x G G} denote the algebraic 
span of all the S^- Then there exist locally normal states on 2IJ(f)o) exhibiting BEC 
for the pure hopping Hamiltonian H = \\A\\% — A. 
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Indeed, consider any Perron-Frobcnius eigenvector u which exists by compact- 
ness, see e.g. Proposition 4.1 of [8]. Fix a > and put 

(5.1) ipaiaH'U'i)aiu2)) -.^ {{e" -iy^ui,U2) +a{ui,u){u,U2) , ui,M2ef)o- 

By transience (cf. (i) of Proposition 4.4 which works for any network G as above), 

(5.1) is meaningful and defines a quasi-free state on 2n([}o). Such a state can be 
interpreted as exhibiting BEC for the presence of a non trivial amount of condensate 
on a lowest energy wave-function u described by the addendum a{ui,u){u,U2). 

As explained before (see also [7]), the careful investigation of the finite volume 
behavior of the KMS states given in (4.8) plays a relevant role because of the 
boundary effects which cannot be avoided in the case of non amenable graphs 
as in our situation. Due to such boundary effects, this analysis becomes very 
delicate. Yet, the emerging situation is surprisingly in accordance with previous 
results in [6, 8]. The main aim of the present section is to investigate, technicalities 
permitting, the infinite volume limits of finite volume Gibbs states arising from 
the grand-canonical ensemble. We start by considering standard finite volume 
approximations of KMS states. We focus the attention to any uniformly bounded 
degree connected network G, equipped with the exhaustion {A„}„gN made of finite 
volume regions. Fix a common root G A„, n G N and consider the sequence 
{i'n}neN of Perron-Frobenius eigenvectors of the Adjacency of the finite regions A„, 
all normalized to 1 at 0. In addition, we also suppose that Aq (or equivalently the 
pure hopping Hamiltonian H = \\Aa\\i — Aa) admits the IDS w.r.t. the exhaustion 
{A„}„gN- We reduce the matter to the C*-subalgebra 2t := 2IT (u„gN^^(A„)) = 
2n(f)o), see Section 5.2.5 of [4]. We denote by I„ := I[£2(a„) the identity on £^(A„). 
Fix a sequence of finite volume chemical potentials {/injneN, with 

f^n< \\Ag\\-\\AaJ, 

converging to /i < 0. The finite volume Gibbs states ujn are given, for such a 
sequence of chemical potential, by 

(5.2) a;„(a^(ui)a(u2)) :== ((e^""^"*" -In)"^ui,M2) , ui,W2et)o. 

The case when /i„ — > /i < is quite standard (cf. [4, 5, 8]) and presents no further 
complications. Thus, we restrict ourselves to the condensation regime /i = 0. In 
correspondence of such a sequence /i„ ^- of chemical potential, we put 

(5.3) Xn ^ \\Ag\\ - ^ln ■ 

The following result takes into account what happens when the graph G is recurrent. 
Proposition 5.2. If G is recurrent then for each x E G, 

\imujn{a^dx)a{Sx)) = +oo . 

n 

Proof. We consider the splitting (4.7), and A„ given in (5.3). It is readily seen 
that A„ > ||Aa„|| and A„ -^ II^gII- By recurrence and the monotone convergence 
theorem, we get 

II X iifc+1 ^ , 1™ „ {RagWSx,Sx) ^ +00 
On the other hand, as j4a„ — > Ac in the strong operator topology, 

{A\jx,Sx) {A^G^x,S,) 

\t' ^ PgII'^+1 ■ 
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Then, again by the monotone convergence theorem and Fatou Lemma, 

+ 00= Hm {RAaWS..S^)^Y. 



K\\Ag\\ 
+ 00 / /Ife 



fc=0 



|A,.||fc+i 



< hm V ^ \2l ^ lim (i?A (A)4,(5^) 



Then we have just shown that the two-point function Ldnia* {6x)a{5x)) diverges in 
the infinite volume hmit. D 



The last result together with Proposition 3.1, tells us that, even when the critical 
density is finite (which certainly happens in presence of hidden spectrum), it is 
impossible to construct locally normal states exhibiting BEC when the graph G 
under consideration is recurrent (compare with the similar results in [8]). In our 
situation, Proposition 5.2 covers the case G'^-'^, with 3 < Q < 7. 

Now we specialize to the situation of main interest here when the network G 
(denoted symbolically by Y) is M^ and G'^''. In this situation, v denotes the 
Perron-Frobenius eigenvector for Ay described in (2.10), (2.11), and w„ the finite 
volume Perron-Frobenius eigenvectors of the finite regions A„. All such Perron- 
Frobenius eigenvectors are normalized at 1 on the fixed common root 0. 

Proposition 5.3. Let w„ be the (unique) Perron-Frobenius eigenvector of Aa„ 
normalized to 1 at the common root G A„. Then lim„ti„(a;) = v{x). In addition 
lim„ ||w„|| ^ +00. 

Proof. We start by considering the networks S„ consisting of G^ , perturbed along 
the subset 5„ :— SnB^, see Fig. 6. The vectors w*^"', also described in [6], are the 




.0- 



1^' 



.0. 



.0. 




Figure 6. The graphs S2 in the cases G^'^ and M^, respectively. 



Perron-Frobenius weights normalized at 1 on the common fixed root. Let r be any 
weak cluster point of {w„}„gN which exists by compactness (cf. Proposition 4.1of 
[8]). It is straightforwardly seen that there exists a collections {fx}x£S of smooth 
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universal functions (on a sufficiently large neighborhood of ||Ay||) such that for 
each X e Ay and n > d{x, 0), 

Here, y{x) € 5* is the unique element such that d{x,S) = d{x,y{x)) (cf. Lemma 
4.2 of [6]), u is either r or v^"\ and, correspondingly, A is \\Ay\\ or ||As^||. Notice 
that, it was shown in [6] that 

limi;*"^(a;) = w(a;) , xeY. 

n 

This first implies that there exists a unique weak cluster point, say r, of the sequence 
{wnjneN- In addition, as ||As^|| '\ ||j4y||, we get for each a; G y, 

v{x) - limi;(")(a;) - lima(|| As„ ||)'^(^'^Vy(x)(PEj|) 

Concerning the limits of the norms we can suppose that the f„ act on £'^{G'^) by 
extending those to on the complement of the balls B^ , respectively. We have by 
Fatou Lemma 

+00 = j|vj| < liminf ||w„|| 

n 

which completes the proof. D 

The precise estimate of the behavior of the i?^-norms \\vn\\ of the Perron-Frobenius 
eigenvectors w„ of Aj^^ as A„ f ^ is needed in the sequel. Unfortunately, even if we 
have a precise estimate of the behavior of k;[^2(y^^J , this does not guarantee the 
analogous estimate (2.12) for ||w„||. However, it is expected that 

(5.4) lim"-"f'f"'^0 

|A„| 

still holds. For the reader convenience, we report a sufficient condition, tested on 
simpler models, which implies (5.4). Suppose that 

Wn(a;) < v{x) , X e 5„ , n e N . 

We automatically get ||u„[£2(5^)|| < ||ti[£2(5^-)|| in i^{Sn) C ^^(C^). Thus, as in 
the proof of Proposition 2.3, we compute, 

II II 2 

<(^^\\u (\\A ihlh2 rrp(s„) \Sn\ ^ 
<(1 + ||^Ab„(||^a„||)||) i^^-n ^U 

because ||z;[£2(5.„)||^|S'„|/|A„| -^ 0, and ||i?AB„ (||^a„||)|| ^ ||i?A^Q (||Ay||)|| as 
Ab„ -^ Aqq in the strong operator topology. The estimate (5.4) implies the weaker 
one 

(5.5) lim ^^ - 

" |A„| 

needed in Proposition 5.7. 

Now we pass to the transient situation by showing a preliminary result which, in 
this case, takes into account the case when the parameter D describing the amount 
of the condensate in (4.9) on the base-point subspace is zero. 
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Proposition 5.4. For each x E Y we have 

lim ((e^" - l)-'5,,5,) = ((e^ - l)-'6,, 5,) . 

Proof. By using the splitting in (4.7), we can reduce the matter to the analysis 
of the Resolvents. Notice that (^A\ Sx,Sx) < (A'^ 6x,Sx)- Thus, by Monotone 
Convergence Theorem, we get 

lim(i?^,„(||A^||)<5.,4) ^limgi^l^ 

+ 00 



J:^-§^ = {RaA\\Ay\\)S.,S.) 



fc=0 



D 



We now recall some notations standardly used in the sequel. Put 

eventually acting on i'^{Sn)- Consider B^ the ball of radius n of G"^, Sn '■— Bn S 
is the support of the perturbation inside A„ (i.e. Sn ^ G'' for G''^'^, and S'„ is 
isomorphic to the finite set consisting of n + 1 points for M.^). Finally, Qn is the 
selfadjoint projection onto i'^(A„) G Cw„. As A„ is the finite volume perturbation 

oiBQ,\\A^c^\\<\\AAj. 

Lemma 5.5. If X e (||AA„|j,r], then 

1 ^r-iiARji 



1-||5„(A)|| - A-PaJI ■ 

Proof. The proof immediately follows from the same computations in Lemma 3.1 
of [6], by taking into account that the resolvent is differentiable and convex. D 

The key-point to investigate the infinite volume limit for the remaining cases, 
that is when the underlying perturbed graph is transient, is the following 

Proposition 5.6. Let ^„ < \\Ay\\ — ||^A„|| (i-e. A„ > mA„||j such that /i„ — t- 

(i.e. A„ — > 1 1 Ay 1 1 j. Then, 

PA„(Mn) - Pc = an+Tn (Q„i?AB„ (A„)-Rs„(A„) (1)-Rab„ (A„)) 

1 

^|A„|(||Ay||-PAj|-M„)' 

where o„ — > 0. In addition, if 

\^ ,0 

\A^\{\\Ay\\ - \\AaJ - fin) 

then phSti-n) -^ Pc- 

Proof. By using the splitting (4.7), the Krein formula (2.6) and finally Lemma 2.4 
of [6], we easily obtain for pA„(Mn) ~ Pc the above formula with a„ given by, 

an =r„ {fiiWAyW - M„)I„ - AaJ) - r^^ (/(Py||I„ - Ay)) 

+r„ (i?^, JA„)) - rA=« (i?^^.,, (II Ay ID) , 

where the last goes to zero by Propositions 2.1 and 2.2. 
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By Lemma 5.5, the second half directly follows as 

r„ {Ra,^ (A„)i?s„(A„)(l)i?A,„ (A„)) 



<||i?^, (PgQ 11)11 V„(i?s„(A„)(l)) 

\S,MAy\\ - \\AbJ) 



<\\RA^i\\AGQ\\)\\ 



\aMAy\\ - \\AaJ\ - ^,r.) - 

D 



Now we specialize the situation to the transient situation for the case when finite 
volume Gibbs states are prepared in order to describe in the infinite volume limit, a 
mean density p > pc. To this end, we chose the sequence of the chemical potential 
by putting 

Even if it is expected that (5.6) implies that the resulting limit density describes a 
mean density p = p^ + c, in order to prove that we should show that 

(5.7) T„ {QnRAn^ (An)il's„(A„)(l)i?A«„ (A„)) ^ . 

The last result (5.7) would follow again by the Krein formula (2.6) for the resolvent. 
Unfortunately, its proof relies by unavoidable boundary effects which seriously affect 
the related formulas and are not easily manageable to take the infinite volume limit. 
However, for the infinite volume mean density p, if (5.7) holds true, we get by 
Proposition 5.6, 

p = limpA„ {p„) ^ pc + lim . ,,|, |, ^-j — ^ r = Pc + O pc ■ 

This simply means that, if (5.7) holds true, the sequence (5.2) of finite volume 
Gibbs states obtained with the choice (5.6) of the chemical potentials, describes 
a fixed mean density p > pc in the infinite volume limit. It is expected that the 
condition (5.6) is incompatible with local normality as it is shown in the following 

Proposition 5.7. Suppose that the estimate (5.5) holds true. If the sequence of 
the finite volume chemical potentials satisfies (5.6), then for the two-point function 
in (5.2), we get 

limu;„(a' ((5o)a((5o)) = +oo . 



Proof. We have 



UJn{a^ iSx)a{Sa;)) > 



\vJ^{\\Ay\\ -||AaJ|-^„) ■ 
The proof easily follows collecting together (5.5) and (5.6). D 

The last result together with Proposition 3.1, would mean that it is impossible 
to construct locally normal states exhibiting BEG with mean density p > p,, in all 
the situations under consideration in the present paper. 
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Finally, we pass for the transient cases to the construction of the states (4.8) 
by the infinite volume limits of finite volume Gibbs states in (5.2) with a suitable 
choice of the sequence of chemical potentials. Also in this situation, we assume that 

(5.8) hm(i?^^JA„)ii's„(A„)(l)i?A^„(A„)g„(5„Q„4) 

= (i?A,(||Ay||)i?s(||^^||)(l)i?^^(||Ay||)<5,,,5,) ,x e G« 
holds true. In order to recover (4.8) by the infinite volume limits of (5.2), we put 

ll«n|P(||Ay|| -||AaJ|-M«) 

Proposition 5.8. Let Y be either M.'^ with 3 < Q < 7, or G'5'« where q > 3, 
q < Q < Q{q) with Q{q) given in (2.8). Suppose that (5.4) holds true, and (5.8) 
is satisfied under the choice (5.9). Then for each x Cz Y , we get for the two-point 
function in (5.2), 

\imuJn{a\Sx)a{dx)) = ujD{a\Sx)a{6,^)) . 

n 

In addition pA^iPn) -^ pc- 

Proof. The second half directly follows from the assumption (5.4), by the second 
half of Proposition 5.6. Indeed, we get, 

\Sn\ ^ l^nlHt^nll' 1 

\AMAy\\ - UaJ - ^in) |A„| ||«„||2(||Ay||-||AAj|-Mn) 

Concerning the first part, as usual we reduce the matter to the Resolvent. As 
A„ -^ W^yW, for each fixed A > ||^y||, we get by monotony 

< (i?A,(A),5„4) - II 112'^"^''^. II, < {RA^AK)QnS.,QnS^) ■ 
By Proposition 5.3, we get after taking the limits in n and in A ^ ll^ylL 

(i?Ay(||^y||)<5:E,4) < hmSUp(i?AA„(-^n)Q«<^a;,Qn4) • 

n 

On the other hand, by using the Krein Formula, we obtain 

{RAA„i\i)QnSx,QnSx) < \RAf.Q {\i)^x, ^x) 
+ {RAB„i^n)Rs,,(\,,)i^)RAB„i^n)Qn^x,QnSx) i 

which, taking into account (5.8), leads to 

liminf (i?AA„(A„)Q„4,Qn<5a;) < (i?A„Q (||^y||)da:, 4) 
+ {Ra^{\\Ay\\)Rs(\IaM^)Ra^{\\Ay\\)5x,Sx) = (i?A,(||Ay||)4,4) . 

Collecting together, we get the assertion. D 

We end by noticing that the case without condensate, but at the critical point 
fi — 0, is directly covered by Proposition 5.4 without any further assumption. 



By transience, it might be enough to check (5.8) only for a fixed point xq in order to obtain 
(5.8) itself for each other point x & Y. 
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